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Abstract 



The coherent and incoherent channels of the neutrinoless muon to electron conversion 
in nuclei, n~ (A,Z) — * e~~ (A, Z)* , are studied throughout the periodic table. The relevant 
nuclear matrix elements are computed by explicitly constructing all possible final nuclear 
states in the context of the quasi-particle RPA. The obtained results are discussed in view 
of the existing at PSI and TRIUMF experimental data for 48 Ti and 20S Pb and compared 
with results obtained by: (i) shell model sum-rule techniques (ii) nuclear matter mapped into 
nuclei via a local density approximation and (in) earlier similar calculations. 
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I. INTRODUCTION 



The neutrinoless muon to electron conversion in the field of a nucleus, 

H~ + (A, Z) - e" + (A, Z)* (1) 

is forbidden in the standard model by lepton flavour conservation and plays an important role 
in the study of flavour changing neutral currents which violate muon and electron numbers 0- 
TTJ. Within the last decade, experiments at PSI ]T|, |3], M and TRIUMF aiming at a search of 



fj, — e conversion electrons have not yet observed such events. These experiments have, however, 
provided us with useful constraints for the violation of muon and electron numbers. The best 
upper limit on the branching ratio 

R T ; e = r(/i~, e -)/r(//- ^) < 4.3 x itr 12 (2) 

has recently been set by SINDRUM II at PSI |I[ by using 48 Ti as target. This value is of the same 
order with the previous limit set at TRIUMF 0, i.e. R T ^ e < 4.6 x 10 -12 . This year |4]], experimental 
data extracted at PSI by using 208 Pb have yielded an upper limit R^ < 4.6 x 1CT 11 . This 
experiment improved by an order of magnitude over the previous upper limit, R^ < 4.9 x 10~ 10 , 
extracted from preliminary experimental data for the same target at TRIUMF B . 



The experimental sensitivity is expected to be further improved by two to three orders of 
magnitude by on going experiments at PSI (to KT 14 ) 0, at TRIUMF (to 10~ 14 ) g § and at INS 
(to 10~ 14 — 10 -16 ) JTTJ . Hopefully, such experiments will not only yield a still better limit, but they 
will detect some (/z~, e~) events which will signal the break down of the muon number conservation 
revealing "new physics" beyond the standard model. For a discussion of lepton flavour violation 
limits in conjunction with theoretical predictions, the reader is referred to the recent survey by 
Depommier and Leroy |§. 

Due to the similarity of electrons and muons, \l — e conversion was originally expected to 
proceed real fast. From a theoretical point of view, the basic background for the (fi~,e~) has 
been set long time ago by Weinberg and Feinberg [|] who assumed that this process is mediated 
by virtual photons (Fig. 1(a)). Non-photonic contributions (Fig. 1(b)- (d)) have been included 
later on in the post gauge theory era (for a recent review on this topic see Ref. ]10[ and for the 
experimental data extracted from various targets see Ref. fll2"| . 



One expects that the Z-exchange diagrams, Fig. 1(b), (c), to be less important than the W-box 
diagrams, Fig. 1(d), even for the incoherent process. The precise value depends, of course, on 
details like quark masses etc. One also expects the W-box to dominate by large factors especially in 
the case of heavy intermediate neutrinos |T3|] . Similar conclusions have been obtained by Marciano 



and Sanda [|Tj]. In the present work we have not included exotic particles like Z' fl5|| , exotic Higgs 
scalars, many Higgs doublets, R-parity violating interactions etc., which may in some models be 
important. We intentionally stayed within the context of the minimal extensions of the standard 
model keeping also in mind that our emphasis here is on the nuclear structure aspects. 

An interesting feature of the (/z~, e~) conversion in nuclei is the possibility of the ground state 
to ground state transitions. The strength of this channel is expected to be enhanced because 
of the coherent contribution of all nucleons of the participating nucleus or at least all protons. 
The rate for such transitions can be expressed in terms of the proton and neutron nuclear form 
factors || [7], |i~6"| . Earlier estimates for the branching ratio R^ e by Weinberg and Feinberg || have 



indicated that, for A > 100 this ratio is approximately constant, while Shanker || found that the 
ratio Rfj, e could be bigger in heavy nuclei. 

The incoherent rate is much harder to calculate. The first such calculations have been per- 



formed only recently |9], DJ in nuclei with closed (sub) shells throughout the periodic table by 
employing shell model sum-rules, i.e. by invoking closure approximation in some form with a 
suitable choice of a mean excitation energy, using a single Slater determinant for the initial state. 
In reality, these calculations give the total rate. The incoherent strength can be estimated by 
subtracting from the total strength the coherent part obtained independently. What, however, 
is needed is the ratio of the coherent rate to the total rate. Shell model results showed that the 
coherent channel dominates the (fi~ , e~) process for light and medium nuclei, but in the region of 
208 Pb, a great part of the rate comes from the inelastic channels. Furthermore, these calculations 
showed that the dependence of the branching ratio on the nuclear mass A and charge Z 
reaches a maximum around A ~ 100 in agreement with the estimates of Ref . fj . 



Recently we have employed, for the coherent and incoherent (n , e ) conversion, another 



approach based on a local density approximation in conjunction with a relativistic Lindhard 
function for the description of the elementary processes: fi~p — > e~p and jjTn — > e^n. The 
incoherent rate in this method was obtained by integrating over the excited states of a local Fermi 
sea. These results have shown that, the coherent contribution is dominant for all nuclei and that 
the branching ratio _R Me presents a maximum in the region of very heavy nuclei, i.e. in the Pb 
region. 



In yet another recent theoretical study of the (fi , e ) conversion [119 , the quasi-particle RPA 



(QRPA) was employed for the construction of the final nuclear states entering the coherent and 
incoherent rate. Such quasi-particle RPA results for A8 Ti have shown that the coherent channel 
dominates. One of the advantages of the QRPA method is that it can be used to estimate the 
mean excitation energy of the nucleus of interest which in turn is useful in checking the results of 
the above mentioned closure approximation which are sensitive to this property. The important 



result |L9[ was that the mean excitation energy E of the nucleus in process (1) is very small, 
E « lMeV, which is appreciably smaller than E « 20MeV used in shell model calculations 
||. The latter value had been chosen from the phenomenology of the charge changing [fi",!/^) 
reaction. This difference is mainly due to the fact that, the coherent elastic channel, possible only 
in the (/i~,e~), is the dominant channel. 

From the above discussion it is clear that, a detailed study of all possible channels of the 
(n~ , e~) conversion for medium and heavy nuclei and in particular for nuclei around 208 Pb, which 
is of current experimental interest, is needed. In the present work, we use the formalism developed 
in the context of quasi-particle RPA [plj (improved in the part of the reduced matrix elements, 
see appendix) to extend our previous results for 48 Ti and we report calculations performed for all 
individual (/i~,e~) conversion channels, in a set of isotopes covering the above region (see below 
Table I) including, of course, 48 Ti and 208 Pb, since, the upper limit on the branching ratio has 
been extracted [Q-|4| from experimental data on these nuclei. We note that, the method using 
a local density approximation []TB[ cannot give us the individual contribution of each accessible 
channel. 

Before embarking on such calculations, we mention that for certain nuclei, in particular those 
with closed shells, like 60 Ni and 208 Pb, a special treatment in QRPA is required in order to 
determine the pairing parameters for protons (gp air ) and neutrons (gp air ). In this work we follow 
the manner used recently in the double beta decay f20| . In sect. II we briefly discuss the method 
used, while the obtained results are presented and discussed in sect. Ill and the conclusions are 
summarized in sect. IV. 



II. BRIEF DESCRIPTION OF THE METHOD 

A. The (jj,~,e~) conversion effective operator 



The effective Hamiltonian operator of the (/i , e ) conversion, which involves both vector and 



axial vector currents, after the usual non-relativistic reduction takes the form [1C 



V v = ^E(3 + fvir:^). fl A = -g A fAY,(i + ""' (3) 

The parameters gv, §a and (3 depend on the assumed mechanism for lepton flavour violation 
0, [TOfl . For the photonic mechanism these parameters take the values gv = 1/6, g A = 0, (3 = 3, 



fv — 1 while, for the non-photonic neutrino mediated mechanism, they are gy = g A = 1/2, 
/3 = 5/6, / y = 1, f A = 1.24 and ^ = /v//a = 1/1-24. 

In Eq. (^|), q represents the momentum transfer to the nucleus. Its magnitude is approximately 
given by 



q = | q | = m M - e 6 - (E f - E gs ) (4) 

where m M is the muon mass, e?, is the muon binding energy and £"/, E gs are the energies of the 
final and ground state of the nucleus, respectively. We should mention that e&, although negligible 
in light nuclei, can become important in heavy elements (see below Table I). 

The matrix elements of the operators of Eq. (|3|) can be obtained via the multipole operators 
Tm^ J given in the appendix (for details see Refs. pi [HI). 



B. Nuclear matrix elements for the coherent rate 



In the case of the coherent (fi , e ) process, 
(0 + — > + ), only the vector component of the (/i~ 



i.e. ground state to ground state transitions 
e~) operator contributes and the coherent rate 



is proportional to |1C 



\<f\n(q) |i,/i>r=& (3 + /W3) 



(5) 



where 



F p ,n(Q 2 ) = I d 3 x Pp , n (x) e- iq - x $ M (x) (6) 

In the last equation p p (x), p n (x) represent the proton, neutron densities normalized to Z and N, 
respectively and $ M (x) is the muon wave function. If we assume that the muon is bound in the Is 
atomic orbit which varies very little inside the nucleus, we can factorize the muon wave function 
out of the integral of Eq. @ and write 

F p (q 2 ) ^< $ ls > ZF z (q 2 ), F n (q 2 ) ^< $ ls > NF N (q 2 ) (7) 
where Fz (-Fjv) the usual proton (neutron) nuclear form factors. 



We should mention that, experimentally the most interesting quantities are the branching ratio 
R^ e and the ratio 77 of the coherent to the total (fj,~,e~) conversion rate (see Eq. flT5| ) below). 
We do not expect the branching ratio to be greatly affected by the approximation of Eq. (|7|), 
especially if we calculate the total /i~ capture rate in the same way. We expect this to be good 
even if for the total muon capture rate we use the Primakoff function which is obtained by 
explicitly using this approximation. The Primakoff function fits the experimental data remarkably 
well throughout the periodic table (even for heavy nuclei). Furthermore, and for similar reasons, 
we expect that the ratio 77 is not going to be drastically affected by this approximation. 

In the above approximation the nuclear dependence of the rate for the coherent process, is 
proportional to the matrix element 



M 2 coh {q 2 ) = Ml^M 2 ) 



1 + 



3-f v /3N F N {q 2 



-1 2 



3 + f v /3 Z F z (q* 



Z 2 Fl(q 2 ) 



(8) 



Thus, the variation of the coherent (fi~, e~) conversion rate through the periodic table can be 
studied by calculating the matrix elements M 2 s ^ gs of Eq. (§) for various A and Z. The nuclear 
form factors involved in M 2 ^„ Q can either be calculated by using various models as shell model 



[p~6| , |22|| , quasi-particle RPA |L9], |23[ etc., or can be obtained directly from experiment whenever 
possible 0, |5|. 



In the context of the quasi-particle RPA with an uncorrelated vacuum as ground state, the 
nuclear form factors are given by (see appendix) 



F T (q 2 



-E 07) (2J + 1) < 3 I Jo(gr) I j >, r = Z, N 



(9) 



where are the occupation probabilities for the proton, neutron single particle states | j 

included in the used model space (j = (n,l,j)). 



> 



We should mention that, in the photonic case (J3 = 3) only the protons of the considered 
nucleus contribute and the right hand side of Eq. (§) becomes Z 2 F^(q 2 ). 



C. Incoherent rate by explicit calculations of the final states 



The incoherent (/i - , e~) conversion rate is evaluated by summing the partial rates for all final 
nuclear states | / > except the ground state. We need calculate the matrix elements for both the 
vector and axial vector operators of Eq. (|3|), i.e. the quantities 

^ = E(— ) 2 /^|</|fia|^>| 2 , \f>*\gs>, a = V,A. (10) 

(q/ is the unit vector in the direction of the momentum transfer qj). 

As we have mentioned in the introduction, for the calculation of Sy and Sa, one can either use 
closure approximation (in which case the state | / >=| gs > is included in Eq. fllOD) or compute 
state-by-state the partial rates involved if one can construct the final states | / > in the context 
of some nuclear model. By using the multipole expansion operators T^ l ^ J (see appendix), the 
matrix elements Sy and Sa are written as 



s a = J2 (^) 2 E l< /-» II f^ J || ^ >| 2 (ii) 

(a = V, A, for the vector (a = 0), axial vector (er = 1) component, respectively). The partial 
matrix element from the initial state + to an excited state | / > in the context of QRPA takes 
the form 



< / 1 1 f^ J 1 1 0+ >= W J X [ X{ MT) U^V}; ] + Y^vVlfg ] (12) 

A,T 

where Vj and represent the probability amplitudes for the single particle states to be occu- 
pied and unoccupied, respectively. They are determined by solving the BCS equations iteratively. 
X and Y represent the forward and backward scattering amplitudes. They are obtained by solv- 
ing the QRPA equations. The index A runs over two particle configurations coupled to a given 
J, namely (ji, j 2 )J for the proton (r = 1) or neutron (r = —1). The quantities W'l = are 
given in the appendix. 

For the total (fi~, e~) rate the relevant matrix elements are obtained by adding the vector and 
axial vector contributions of the coherent and incoherent rate i.e. 



M 2 0t = S V + 3S A + So (13) 
where Sq is associated with the ground state to ground state transition 



for the vector component a = 0, and for the axial vector component a — 1. For + nuclei only 
the vector term contributes. 



III. RESULTS AND DISCUSSION 



Using the method outlined above, in the present work we have calculated the matrix elements 
needed for both the coherent and incoherent (fi~,e~) rate, for the nuclei AS Ti, m Ni, 72 Ge, 112 Cd, 
162 Yb and 208 Pb. The specific parameters used and a brief description of the model spaces employed 
can be read from Table I. For all nuclei considered we have employed the same model space for 
protons and neutrons. For the harmonic oscillator parameter b in the region of heavy nuclei we 



have employed the improved expressions of Ref. [26 



In the BCS description of the uncorrelated ground state, for each nuclear isotope the single 
particle energies have been calculated from a Coulomb corrected Wood-Saxon potential with spin- 
orbit coupling. The G-matrix elements of the realistic Bonn one-boson exchange potential P?| 



have been employed. The values of pairing parameters, g v vaiT and gp air , renormalizing the proton 
and neutron pairing channels in the G-matrix have been deduced by comparing the quasi-particle 
energies with experimental pairing gaps as is described in Refs. |]28| , |29| . For the special cases of 
2gAz, which is a proton closed-shell nucleus and ]^Pb, which is a doubly closed-shell nucleus, the 
pairing parameters have been deduced from the neighboring nuclei ^Fe and fl 8 Po, respectively, 



in analogy with the procedure followed in the study of the nuclear double beta decay in the double 
closed shell nucleus 20^° P0|] - The resulting pairing parameters gp air , g^ aiT for each nucleus are 
shown in Table I. 



A. The coherent process 



It is obvious from Eq. (H) that, for the coherent process, i.e. gs —>■ gs transitions, we need the 
proton and neutron nuclear form factors, F z (q 2 ) and F N (q 2 ), respectively. The results obtained 
by using as ground state the uncorrelated RPA vacuum, are listed in Table II for the following 

two cases: 

(i) By neglecting the muon binding energy e& in Eq. (f|) (as in Refs. 0, Then, the elastic 
momentum transfer is the same for all nuclei, i.e. q ~ m M ps .535/m~ 1 . Such results are indicated 
as QRPA(i). 

(ii) By taking into account in Eq. (f|). Then, the elastic momentum transfer is q ~ m M — e& 
and varies from q ps .529fm~ 1 , for 48 Ti where e?, ps 1.3MeV, to q ~ .482/m -1 , for 208 Pb where 
€b ps 10.5MeV (see Table II, results indicated as QRPA(ii)). 

In Table II we also present the shell model results of Ref. 0, obtained with q = .535 fm~ x 
throughout the periodic table, i.e. as in case (i) above. We see that, QRPA(i) and shell model 
methods, give about the same results. However, the form factors of QRPA(ii) for heavy nuclei 
differ appreciably from those of both QRPA(i) and shell model. For 208 Pb, for example, the 
QRPA(ii) form factors are about 30% larger than the corresponding QRPA(i) and shell model. 
This happens because the inclusion of results in a smaller momentum transfer to the nucleus 
and, consequently, in an increase of the form factors. The larger the value of e&, lead region, the 
bigger the difference between form factors QRPA(i) and QRPA(ii). 

In Table II we also show the experimental form factors obtained from electron scattering data 
24] , |25| at momentum transfer q = m M — e&. We see that, when using the right form factors, 



i.e. taking into account the binding energy e& (QRPA(ii) case), the form factors calculated in the 
present work, are in good agreement with the experimental ones. The deviation is less than 5% 
with the possible exception of 112 Cd and 208 Pb where it is about 10%. 

The variation of the coherent nuclear matrix elements M 2 oh with respect to A and Z is shown 
in Fig. 2(a), for the photonic mechanism, and Fig. 2(b), for the non-photonic one. From these 
Figures we see that, by taking into account the muon binding energy e?,, QRPA(ii), all matrix 
elements increase continuously up to the lead region where they become about a factor of two 
bigger than the corresponding QRPA(i) and shell model values. This implies that the coherent 
rate becomes bigger for heavy nuclei, Pb region, which makes such nuclei attractive from an 
experimental point of view pi [|] provided, of course, that they also satisfy other additional criteria, 
e.g. the minimization of the reaction background etc. || l^j. The , e~) conversion electrons of 



a given target are expected to show a pronounced peak around E e = m M — e b , which for the lead 
region is E e ~ 95MeV. One prefers this peak to be as far as possible above the reaction induced 
background. 

We should recall that, in the present work and in shell model method of Ref. 0, the factor- 
ization approximation Eq. ([?]) was used. The exact expression, Eq. (El), was used in Ref. 



and yielded matrix elements which for heavy nuclei are bigger than the approximate ones. This, 
however, as we have extensively seen in sect. II. B, only slightly affects the branching ratio 
and the ratio rj of the coherent rate to the total rate, which in our case are the most important 
quantities. We also mention that, shell model results for the total muon capture rate pHI], ob- 



tained by using the exact muon wave function, differ by only 5.7%, in the case of 60 Ni, and by 
7.0%, in the case of 208 Pb, from those obtained by using the approximation of Eq. (|7j). Detailed 
QRPA calculations, which do not invoke this factorization approximation, are under way and will 
be published elsewhere. 

B. Incoherent process 



As we have stated in sect. II. C, the incoherent process in the present work is investigated 
by calculating state-by-state the contributions of all the excited states of the nucleus in question 
which are included in the model space described in Table I. 

For the photonic mechanism the nuclear matrix elements obtained for all positive and negative 
parity states up to 6~, 6 + are shown in Table III. For this mechanism only the vector component, 
Sy gives non-zero contribution (Mf nc = Sy)- For the non-photonic mechanism we have non zero 
contributions from both the vector and axial vector components, Sy and Sa, and the results are 
shown in Table IV (Mf nc = S v + 3S A ). 

From Tables III and IV we see that, the main contribution to the incoherent rate comes from 
the low-lying excited states. High-lying excited states contribute negligibly. This means that for 
a given A, a nuclear isotope with many low-lying states in its spectrum is characterized by big 
incoherent matrix elements. 

For the doubly closed shell nucleus 208 Pb, which is of current experimental interest 0, 
the incoherent matrix elements are smaller than expected. A plausible explanation is that the 
spectrum of this nucleus presents a big gap (minimum energy needed to excite the first excited 
state) and only few excited states lie below ~ 5MeV. 

In general, the incoherent matrix elements do not show clear A and Z dependence. Their 
magnitude depends on the spectrum of the individual nuclear isotope. 

In obtaining the results of Tables III and IV for 1~ states, we removed the spurious center of 
mass contributions by explicitly calculating the spurious state \S > and removing its admixtures 
from the incoherent and total rate. We have also calculated the overlaps < l~,m\S > (where m 
counts the 1~ excited states) and found that most of the spuriousity lies in the lowest 1~ state 
being 88% for m Ti, 63% for 60 A^, 62% for 72 Ge, 57% for 112 Cd, 87% for W2 Yb and 77% for 208 Pb 
3D[. We should stress that, for all nuclei studied, the spurious center of mass contribution is less 



than 30% of the incoherent matrix elements, i.e., 1.0-1.5% of the total (/i , e ) conversion rate. 

An additional point we should note is the effect of the ground state correlations on the e~) 
matrix elements. In QRPA this can be easily estimated by using a correlated quasi-particle RPA 
vacuum instead of the uncorrelated one [|32]-[|35|. In the present work we have not performed 
additional calculations with a correlated RPA vacuum. It is known, however, that the matrix 
elements for 48 Ti obtained this way [[19] are reduced by ~ 30%. The ground state correlations 
tend to decrease the strengths of all (/i - , e~) conversion channels, but do not affect the parameter 
rj (see sect. III.C). 

C. Comparison of coherent and incoherent processes 



As we have seen above, the gs — > gs channel is the most important one. Therefore, a useful 
quantity for the (fi~,e~) conversion is the fraction of the coherent matrix elements M 2 oh divided 
by the total one M 2 ot , i.e. the ratio 



V = MlJMl t 



(15) 



In earlier calculations 77 was estimated H to be a decreasing function of A with a value of 
77 ~ 83 % in Cu region. By using, however, the most appropriate QRPA(ii) results, we find 
that indeed the coherent channel dominates throughout the periodic table (see Table V, for the 
photonic and Table VI, for the non-photonic mechanism). In fact we see that the values of 77 



obtained in the present calculations are a bit bigger than those of Ref. |18| obtained with a local 
density approximation and a lot bigger than those of Ref. ||. We should stress, however, that the 
exaggeration of the incoherent channels in shell model calculations of Ref. |§ is not a shortcoming 
of the method itself but the result of ignoring the muon binding energy e& in calculating the nuclear 
form factors. In fact, repeating the calculations of Ref. || and taking into account the effect of 
e& on the form factors of the coherent process as well as on the mean excitation energy entering 
the total rate, we find a value of 77 > 75%. 

IV. CONCLUSIONS 

In the present work we have studied in detail the dependence of (/j~,e~) conversion matrix 
elements on the nuclear parameters A and Z. Our nuclear matrix elements were obtained in the 
context of Quasi-particle Random Phase Approximation (QRPA) which permits a relatively simple 
construction of all needed final states. So, there was no need to invoke closure approximation. The 
results obtained cover six nuclear systems from 48 Ti to 208 Pb, which are of experimental interest. 
The most important conclusions stemming out of our detailed study are: 

i) The coherent mode dominates throughout the periodic table but it is more pronounced in 
the heavy nucleus 208 Pb which is currently used at PSI in the SINDRUM II experiment. 

ii) The coherent and total rates as well as the ratio 77 (coherent to total) tend to increase as a 
function of the mass number A up to the Pb region. This encourages the use of heavier nuclear 
targets to look for lepton flavour violation. 

iii) In evaluating the nuclear matrix elements the muon binding energy should not be ignored 
especially for heavy nuclear elements. 

iv) The great part of the incoherent rate comes from the low lying excitations. 

The results obtained in the present work are in good agreement with those obtained in the 
framework of the local density approximation |Tj| as well as those of shell model calculations 



provided that all calculations take into account the muon binding energy e^. 
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APPENDIX 

A). The multipole expansion operators T^' a ^ J , resulting from Qy and ft a of Eq. ([|), are 
written as 



p(l,0)J 
■ M 



g v 5 u V^J2( 3 + ^ 3i )ji(qri)Y^{ri) 



(16) 



i=i 



for fly, the spin independent component, and 



i=i 



M 



(17) 



for Qa, the spin dependent component. ji{qr) are the spherical Bessel functions. 

The quantities W( = Wf^ of Eq. ( |T2| ) contain the reduced matrix elements of the operators 
T J between the single particle proton or neutron states ji and j 2 as 



= (C + T0) 



<3i \\T J \\j2> 
2J+1 



18) 



(( = 3, for fly and ( = 1/1.24, for Qa)- The reduced matrix elements < j'i||T J ||j2 > are given in 
Ref. [l!J. The relevant radial matrix elements < niii|j/(gr)|ri2/2 > ; for harmonic oscillator basis 
often used, can be written in the elegant way 



< nih\ji(qr)\n 2 l 2 > = e"* E ^X K+l/ \ X = {qbf/A 



K=0 



where 



(19) 



ft; 



= m + n 2 + m, m = (h + l 2 — l)/2 
The coefficients e K (nili,n 2 l 2 ,l), in general simple numbers, are given by 



nni\n 2 \ 



4r> 1 + z 1 + f)i> 2 + / 2 + f) 



i ni ?12 
2 



E E n!A Kl (n 1 / 1 )A K2 (n 2 / 2 )A K (n/) (20) 



where the A K (nZ) are defined in Ref. [[[£]], n — Ki + K 2 + m and 



0. 



k — m — n 2 < 



k — m — n 2 , k — m — n 2 > ' 



0, k — m — Ki < 

k — in — Kx, k — m — Ki > 



The advantage of Eq. (|19|) is that, it permits the calculation of e K , which are independent of 
the momentum q, once and for the whole model space used. Afterwards, the relevant reduced 
matrix elements are easily obtained for every value of the momentum transfer q. 



B). In the context of the quasi-particle RPA, the point-proton (-neutron) nuclear form factors 
of Eq. (|5|) can be cast in the compact form 



F T (q 2 ) = -e- W 2 / 4 ff 61 (qbf\ r = Z, N (21) 



r A=0 



where b is the harmonic oscillator parameter, N space represents the maximum harmonic oscillator 
quanta included in the model space used (see Table I), and 6\ the coefficients 



v^F „ , A 2 (2j + l)n!C, 



x-i 
ni 



where {vj^j are the occupation probabilities for the proton (neutron) single particle j-levels. The 
coefficients C™[ are given in Ref. p2| . 
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TABLE I. Renormalization constants for proton (gp air ) and neutron (<7p air ) pairing interactions 
determined from the experimental proton (A? ep ) and neutron (A^ xp ) pairing gaps. 



\ n p Qii c 

in ucieus 


V^OIlllg, UI cLLlUIl OpcLCc 


°ho\J m ) 


A exp f ]\/Tp\f\ 


A expf if„T/\ 


n p 

iJpair 


n 

Upair 


48 TV 

22 1 ? 26 


io levels ^no coiej 


1 Q9 

i.yz 


l.oyO 


1 .0041 


1 .UoZ 


l.UUZ 


l° 8 Ni 32 


16 levels (no core) 


2.02 


1.718 a 


1.395 a 


1.033 


0.901 


32 (j-e4 


16 levels (no core) 


2.07 


1.611 


1.835 


0.924 


0.995 


4g OaQ4 


16 levels (core 2oC a 2o) 


2.21 


1.506 


1.331 


1.099 


0.950 


162 v/, 
70 r °92 


23 levels (core 20^020) 


2.32 


1.170 


1.104 


0.894 


0.951 


208 pi. 
82 - r °126 


18 levels (core Iq S71 50 ) 


2.40 


0.807 a 


0.611 a 


0.861 


1.042 



a For the closed shell nuclei the parameters g v vair and g^ air have been borrowed from the (N ± 
2,Z=)=2) nuclei i.e. the experimental gaps (columns 4 and 5) for^Ni^ andf^Pbue, are those of 
26-^ e 34 an d 84 8 -P°i24; respectively. 



TABLE II. Nuclear form factors for protons (F%) and neutrons (-Fat) calculated in the context 
of the shell model |J and quasi-particle RPA cases: QRPA(i) and QRPA(ii) (see text). For 
comparison the experimental form factors |24| ^] are also shown. 



Nucleus 


Shell Model 




QRPA(i) 


QRPA(ii) 




Exper. 


(A,Z) 


ho (fm^ 1 ) 


F z 


F N 


F z 


F N 


e b (MeV) 


F z 


F N 


pexp 


22^26 


1.906 


.543 


.528 


.528 


.506 


1.250 


.537 


.514 


.532 


t%Ni 32 


1.979 


.489 


.478 


.489 


.476 


1.950 


.503 


.490 


.494 


llGe^o 


2.040 


.470 


.448 


.456 


.435 


2.150 


.472 


.451 


.443 


112^J 
48 L-«64 


2.202 


.356 


.318 


.349 


.312 


4.890 


.388 


.352 


.353 


70 1 °92 


2.335 


.261 


.208 


.252 


.218 


8.445 


.314 


.280 


.305 


208 pi, 
82 - ro 126 


2.434 


.194 


.139 


.207 


.151 


10.475 


.271 


.214 


.242 



TABLE III. Incoherent fi — e conversion matrix elements (Mf nc ) for the photonic mechanism. 



r «Ti SSJV* gGe ii 2 Cd i^yfe 1S 8 P& 



0+ 1.946 1.160 2.552 2.088 4.305 2.512 

2+ 0.242 0.738 1.396 2.669 6.384 2.342 

4+ 0.004 0.005 0.015 0.021 0.063 0.056 

6+ 6 10~ 6 6 10~ 6 1 10" 5 6 10~ 5 2 10~ 4 3 10~ 4 

1- 3.711 4.215 5.066 5.282 4.824 4.533 

3~ 0.037 0.081 0.152 0.249 0.542 0.476 

5" 2 10~ 4 2 10" 4 5 10~ 4 0.001 0.005 0.005 

S v 5.940 6.199 9.181 10.309 16.123 9.924 



ML 5.940 6.199 9.181 10.309 16.123 9.924 



TABLE IV. Non-photonic mechanism. Incoherent \x — e conversion matrix elements (M? nc ) for 
the vector (Sy) and axial vector (Sa) component of the (/i~,e~) operator. 



r f,Ti ™Ni llGe \?Cd % 2 Yb ™Pb 
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2 10~ 4 


4 10~ 4 
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6.796 


5.537 
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0.328 


0.684 


0.572 


5" 


1 10" 4 


2 10" 4 


5 10" 4 


0.002 


0.006 


0.006 


S v 


6.679 


8.708 


12.095 


15.477 


26.280 


14.061 


1+ 


0.265 


1.114 


0.795 


0.943 


1.599 


0.951 
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0.041 
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0.362 


0.644 


0.466 
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0.221 


0.244 


0.422 


0.693 


0.635 


4+ 


2 10~ 4 


0.001 
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0.005 


0.013 


0.015 


5+ 


2 10~ 4 


0.002 
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0.013 


0.021 


6+ 


1 10- 7 


i io- 6 


1 io- 7 


1 10~ 5 


4 10~ 5 


9 10~ 5 


0" 
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1.498 


1.394 
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1.539 


\- 
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1.137 


1.170 


3" 


0.010 


0.015 


0.024 
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9 10- 6 


7 10- 5 


1 io- 4 


3 10~ 4 


0.001 


0.001 


3S A 


2.405 


5.009 


4.637 


6.115 


7.305 


6.226 



ML 9.084 13.717 16.732 21.592 33.585 20.287 



TABLE V. Coherent and total fi — e conversion rate matrix elements QRPA(ii) (see text) for the 
photonic mechanism in the neutrino mediated process. For comparison, we also show the ratio 77 
of Eq. ( |15"D given by shell model obtained by ignoring the muon binding energy e^. The results 
of QRPA(i) are similar to those of the shell model. 



AT 1 

Mucleus 


QRPA(n) 


TIT 1 ■ TH1 j_ 

Matrix Elements 




r]% 


(A rv\ 

(A z) 
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gs^gs 


M tot 


QRPA(n) 


r~i n 11 n r 1 1 

Shell Model 


48 T 1 ' 
22 1 ? 26 


139.6 


145.5 


95.9 






198.7 


204.9 


96.9 


64.9 


3 |Ge4o 


227.8 


237.0 


96.1 


59.7 


4g Ott64 


346.7 


357.0 


97.1 




70 r °92 


484.3 


500.4 


96.8 


36.9 


208 pi. 
82 r °126 


494.7 


504.6 


98.0 


25.5 



TABLE VI. The same as in Table V for a non-photonic mechanism (/3 = 5/6) in the neutrino 
mediated process. For comparison we have added the results for rj obtained with Local Density 
Approximation (L.D.A.) [|18| and shell model ||. 



Nucleus 


QRPA(ii) 


Matrix Elements 




v% 




(A,Z) 


M 2 


ML 


QRPA(ii) 


Shell Model 


L.D.A. 


22^26 


375.2 


384.3 


97.6 




91.0 




527.4 


541.1 


97.5 


74.9 






639.5 


656.2 


97.4 


70.3 




\Y2fi J 
48 <~^«64 


983.3 


1004.9 


97.8 




93.0 


70 r °92 


1341.2 


1374.8 


97.6 


40.1 




208 pi. 
82 - r °126 


1405.2 


1425.5 


98.5 


28.2 


94.0 



Figure Captions 



FIGURE 1. Typical diagrams entering the (/i , e ) conversion: The photonic (a), Z-exchange 
(b),(c) and box (d) diagrams. Only the specific mechanism involving intermediate neutrinos is 



exhibited here. v e = J2j UejVj, 



J2j UujVj, where Uj are the neutrino mass eigenstates and U, 



<-j ■ 



U^j are the charge lepton current mixing matrix elements. Other mechanisms can also contribute 
(SUSY, Z', Higgs etc., see Ref. M). 



FIGURE 2. Variation of the coherent (yU~, e~) conversion matrix elements M^ oh for specific mass 
A and charge Z (see text) for the photonic mechanism (a) and the non-photonic mechanism (b). 
In QRPA(i) the muon binding energy e& was neglected, but it was included in QRPA(ii). We see 
that eb strongly affects the matrix elements for heavy nuclei. For comparison the results of Ref. 
(shell model results) are also shown. For photonic and non-photonic diagrams the coherent 
rate increases up to Pb region where it starts to decrease. 
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